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Electron-hole correlation in quantum-dot quantum wells 
(QDQW's) is investigated by incorporating Coulomb and ex- 
change interactions into an empirical tight-binding model. 
Sufficient electron and hole single-particle states close to the 
band edge are included in the configuration to achieve con- 
vergence of the first spin-singlet and triplet excitonic energies 
within a few meV. Coulomb shifts of about f 00 meV and ex- 
change splittings of about f meV are found for CdS/HgS/CdS 
QDQW's (4.7 nm CdS core diameter, 0.3 nm HgS well width 
and 0.3 nm to 1.5 nm CdS clad thickness) which have been 
characterized experimentally by Weller and co-workers [ D. 
Schooss, A. Mews, A. Eychmuller, H. Weller, Phys. Rev. B, 
49, 17072 (1994)]. The optical excitonic gaps calculated for 
those QDQW's are in good agreement with the experiment. 
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Semiconductor nanostructures, from quantum wells, 
quantum wires to quantum dots 1,2 , have been ex- 
tensively investigated due to their remarkable appli- 
cations, for example, as fast and low noise elec- 
tronic devices and tunable optoelectronic elements. 
Recently, a class of new and promising hetero- 
quantum dots, termed quantum-dot quantum wells 
(QDQW's), have been successfully synthesized in wa- 
ter, for example, CdS/HgS/CdS 3-5 , CdTe/HgTe/CdTe 6 
and ZnS/CdS/ZnS 7 . These QDQW's have internal 
nanohctcrostructures with a quantum-well region con- 
tained inside the quantum dot. Meanwhile, high- 
resolution transmission electron microscopy images 5 ' 6 
have shown that CdS/HgS/CdS and CdTe/HgTe/CdTe 
QDQWs are not spherical, but are preferentially trun- 
cated tctrahedral particles. However, spherical shell par- 
ticles are commonly considered to explain experimental 
results 3 ' 4 ' 8 " 12 . 

Numerical calculations 3 ' 4 ' 8 ' 9 on QDQW's have been 
based generally on the one-band effective-mass approx- 
imation as well as the parabolic approximation for the 
conduction and valence bands. These theoretical stud- 
ies have demonstrated the remarkable effect of the in- 



ternal well on single-particle electron and hole energies 
and pair overlaps 3 ' 4,8 and determined the ground-state 
energy of an uncorrelated electron-hole pair 3 ' 4 . More- 
over, Bryant 9 determined the contribution of pair corre- 
lation to the electronic structure of QDQW's nanosys- 
tems. Very recently, Jaskolski and Bryant 10 have devel- 
oped a multi-band theory to determine electron, hole and 
exciton states for QDQW's. Actually, an atomic model 
is essential for QDQW's since the internal wells are no 
more than a few monolayers thick 11,12 . Hence, in this 
paper, we incorporate for the first time Coulomb and ex- 
change interactions into an empirical tight-binding model 
to describe QDQW nanocrystals. In detail, we shall in- 
vestigate electron-hole interactions and optical excitonic 
gaps of QDQW's, and report our numerical results by 
considering three CdS/HgS/CdS QDQW's which were 
characterized experimentally by Weller and coworkers 3 . 
Comparison with the experiment shows that our tight- 
binding theory provides a good description for QDQW 
nanocrystals. 

First, we use the empirical tight-binding (ETB) 
method 13-16 to perform numerical calculations of both 
electron and hole single-particle energies (i.e, E e and Eh) 
and eigenstates (i.e., | $ e ) and | <&/j}) for QDQW's. Our 
theory can be used to model single or coupled nanocrys- 
tal systems with spherical, hemi-spherical, tetrahedral or 
pyramidal geometry. In this paper, we focus on single 
spherical QDQW's. Also we assume that atoms in these 
nanoparticles occupy the sites of a regular fee lattice. As 
developed by Vogl, Hjalmarson and Dow 17 , each atom 
has its outer valence s orbital and three outer p orbitals 
plus a fictitious excited s* orbital which is included to 
mimic the effects of higher lying states. Only on-site 
and nearest-neighbor couplings between orbitals are in- 
cluded in our sp 3 s* ETB theory. Since the spin-orbital 
coupling is not known for HgS, we shall not consider it 
for our numerical results reported here. The empirical 
single-particle Hamiltonians are determined by adjusting 
the matrix elements to reproduce known band gaps and 
effective masses of the bulk band structures. In this pa- 
per, our tight-binding parameters for CdS and HgS are 
the same as those of Bryant and Jaskolski 11,12 . Finally, 
the electron and hole single-particle eigenstates and en- 



1 



ergies close to the band edges are found by diagonalizing 
the single-particle Hamiltonian with an iterative eigen- 
value solver. 
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FIG.l: Electron and hole single-particle energy spectra, E e 
and E h , for (4.7, 0.3, 0.3), (4.7, 0.3, 1) and (4.7, 0.3, 1.5) 
CdS/HgS/CdS QDQW's. The level degeneracy is shown in 
the brackets and the approximate spherical symmetry of each 
state is indicated. Both the electron energy, which increases 
upward, and the hole energy, which decreases downward, are 
referred to the top of the valence band of CdS. 

We describe the effective Hamiltonian of an electron- 
hole pair by combining a two-particle term, H e h, which 
includes the Coulomb and exchange interactions, with a 
single-particle term, H S i„ g i e , which contains the kinetic 
and potential energies of the electron and hole 16 . The 
Coulomb and exchange interactions are screened by a 
dielectric function e(| r r \,R), where | r — r and 
R are the separation between electron and hole and the 
quantum dot radius, respectively. Our electron-hole basis 
set {| ^>i}} is taken by multiplying the spatial part | § e h) 
(namely, the product of electron and hole single-particle 
eigenstates | $ e ) and | <&/j)) and their spin states | (j> sp i n ) 
(namely, either the singlet component or one of the triplet 
components of the electron-hole spin state). Then, the 
single-particle Hamiltonian can be written in terms of the 
electron-hole basis set as 



where E { = Ef> - E h v is the energy difference between 

corresponding electron and hole energies, Ei l) and £?« of 
the single-particle Hamiltonian, of the zth electron-hole 
pair. Meanwhile, the electron-hole interaction Hamilto- 
nian in the electron-hole basis set is given by 



*i><*i 



(1) 



(2) 



where J and K describe the Coulomb and exchange inter- 
actions (the same as those of Lee et al. 16 ). More details 
about the formulas mentioned in this part are given in 
the work of Lee et al. 16 . 

TABLE I: On-site unscreened Coulomb and exchange inte- 
grals, uJcoui and uj exc h for the sp 3 s* basis set in units of eV 
for Cd, Hg and S. Integrals for the sp 3 orbitals are calculated 
based on the hybridized orbitals along the bonding directions 
defined by Leung and Whaley 14 . 
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TABLE II: Coulomb shift E cou i, exchange splitting E exc h and 
the first spin-singlet and triplet excitonic energies e[^ and 
e[ 3) for a (4.7, 0.3, 0.3) CdS/HgS/CdS QDQW. iV e and N h 
denote the number of electron and hole single-particle states 
included in the basis set. 
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1 3 


2079.62 2079.11 


94.01 


0.51 


4 3 


2073.23 2072.46 


100.66 


0.77 


4 6 


2067.28 2066.22 


106.90 


1.06 


4 9 


2067.19 2066.12 


107.00 


1.07 



As introduced by Leung and Whaley 14 and Lee et al. 16 , 
the Coulomb and exchange interaction matrix elements 
are expressed in terms of the Coulomb and exchange inte- 
grals, uJcoui and uj exc h, of our ETB orbitals. Table I lists 
the unscreened on-site Coulomb and exchange integrals 
for the sp 3 s* basis set for Cd, S and Hg calculated by 
using a Monte Carlo method with importance sampling 
for the radial integrations 16 . Regarding off-site Coulomb 
integrals, we estimate them using the Ohno formula mod- 
ified by Leung and Whaley 14 . It is known that off-site 
exchange integrals decrease quickly as the distance be- 
tween atom sites increases due to the localization and 
orthogonality of orbitals 16 . Here, we use the off-site ex- 
change integrals of Leung, Pokrant and Whaley 15 . 
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The high-frequency dielectric constant 3 for the differ- 
ent regions of CdS/HgS/CdS QDQW's are e cd s = 5.5, 
e H g s — 11-36, and ejf 2 o = 1-78. In our numerical cal- 
culations, we use an average, effective high-frequency di- 
electric constant, e ave = 6. 
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FIG. 2: Spin-singlet excitonic energy spectrum for the first 
few excitonic states (j, the state serial number), where A, • 
and □ arc for (4.7, 0.3, 0.3), (4.7, 0.3, 1) and (4.7, 0.3, 1.5) 
CdS/HgS/CdS QDQW's, respectively. The primary electron- 
hole pair states making up each of the triply degenerate exci- 
tonic levels are indicated. 
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FIG. 3: Occupation probability P(K) of the Kih electron- 
hole pair (n e ,rih) in the first (solid line), 5th (dotted line), 8th 
(dashed line) and 22nd (dot-dashed line) spin-singlet excitonic 
states for the (4.7, 0.3, 0.3) CdS/HgS/CdS QDQW. P(K) is 
broadened by a Gaussian to enhance visualization. 

Because of the total spin of the electron-hole pair, we 
have two Hamiltonians: one for a spin singlet includ- 
ing both the Coulomb and exchange interactions, and 
another for a spin triplet having only the Coulomb in- 
teraction. Therefore, we diagonalize both Hamiltonians 
in the electron-hole basis separately and obtain a set 



of spin-singlet and triplet excitonic states. For clarity, 
we use the parameter set (D core , L we u, L c i a( i) to de- 
note a CdS/HgS/CdS QDQW, where D corei L weU and 
Ldad are the CdS core diameter, HgS well width and 
CdS clad thickness, respectively, in nm. Also we men- 
tion three definitions of Lee et aJ. 16 : (i) optical exci- 
tonic gap E° pt is the lowest spin-singlet excitonic energy 

E[^] (ii) the difference between the single-particle energy 
gap E s g mgle and the lowest spin-triplet excitonic energy 

(3) 

E\ is defined as the Coulomb shift, namely, E cou \ = 
^single _ £>(3) . ^--j ^ ne difr crcncc between the lowest 
spin-singlet and triplet excitonic energies is defined as 
the exchange splitting, namely, E exc h — e[^ — e[ 3 \ 

For the band gaps of HgS and CdS, we use E g ,H g s = 0.2 
eV and E g .cds = 2.5 eV, respectively, with the CdS con- 
duction band edge 1.45 eV above the conduction-band 
edge of HgS 10 . Fig.l shows the electron and hole single- 
particle energy spectra of (4.7, 0.3, 0.3), (4.7, 0.3, 1) and 
(4.7, 0.3, 1.5) CdS/HgS/CdS QDQW's, which have been 
characterized experimentally by Weller and coworkers 3 . 
The calculated single-particle energy gaps are 2.173 eV, 
2.099 eV and 2.085 eV, respectively. The degeneracy is 
also shown in the brackets. The lowest electron and hole 
states can be described approximately by the spherical 
symmetries indicated in Fig.l, i.e., 1S-, IP- and ID-like 
electron states, mixed lPo + lP2-like hole states, and 
a lSi-like hole state. Here, S, P, D indicate the spa- 
tial angular momentum and the subscript for holes indi- 
cates the total angular momentum sum of the spatial and 
atomic orbital angular momenta. As the clad thickness 
increases, the level ordering of the electron states does 
not change, but that of the last two hole states switches 
(as shown in Fig.l). The ordering of hole levels switches 
because excited P hole levels are the states most sensitive 
to the potential far from the dot center. 

In this work, we focus mainly on the lowest singlet 
and triplet excitonic energies. In our electron-hole ba- 
sis set, we include the 4 electron (n e = 0,1,2,3) and 9 
hole (n/j = 0, 1, ...,8) single-particle states closest to the 
band edge, where n e (labelling initially from the ground 
electron state) and nh (labelling initially from one of the 
triply degenerate ground hole states) are the indices for 
electron and hole single-particle states. For example, Ta- 
ble II lists the Coulomb shift, exchange splitting and the 
lowest spin-singlet and triplet excitonic energies for the 
(4.7, 0.3, 0.3) QDQW by considering different numbers of 
electron and hole single-particle states close to the band 
edge. The cases shown are chosen to ensure that all states 
of a given degenerate level are included. The inclusion of 
more electron-hole configurations leads to an increase of 
the Coulomb shift and exchange splitting. The excitonic 
energies are converged within a few meV when N e = 4 
and N h = 9. 

Fig. 2 shows the singlet excitonic energy spectrum for 
the lowest few excitonic states in the three QDQW's. 
For the electron-hole basis, N e = 4 and Nh = 9, 
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the first three triply degenerate excitonic energy levels 
come mainly from the contribution of electron-hole pairs 
formed with the ground-state electron, as indicated in 
Fig. 2. Other higher excitonic energy levels are due to 
electron-hole pairs made from the excited-state electron. 
This is clearly shown in Fig. 3, which presents the occu- 
pation probability P(K) of the Kth electron-hole pair 
(n e ,nh), K = (n e ,nh) = fih + 1 + 9n e , for one exciton 
state from each of the three lowest triply degenerate ex- 
citon levels (the first, 5th and 8th exciton states in Fig. 2) 
and a higher exciton level (the 22nd exciton state). Fig. 3 
shows that correlation does not strongly mix pair states 
of different energy in forming the low excitonic states. 
This explains why convergence is achieved with only a 
few basis states. Our previous calculations of correla- 
tion effects in QDQWs 9 showed that correlation effects 
are weaker in QDQWs than in quantum dots. The extra 
local confinement of the electron and hole to the quan- 
tum well inside the quantum dot increases the splitting 
between the lowest single-particle states and suppresses 
Coulomb mixing of these states. 

TABLE III: The Coulomb shift E cou i, exchange splitting 
E exc h, optical excitonic gap E° pt , and 1S-1S transition en- 
ergy Eis-is calculated for (4.7, 0.3, L c i ad = 0.3,1,1.5 ) 
CdS/HgS/CdS QDQWs. Ef v is the band gap measured by 
Weller's group 3 . EZ"r and EfgHg are the Coulomb shift 
and IS — IS transition energy, respectively, calculated with 
the effective mass approximation. 
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meV to 0.8 mcV when L c i a d varies from 0.3 nm to 1.5 
nm. 

Including a finite barrier to represent the water so- 
lution and the Coulomb interaction between electron 
and hole, Weller and coworkers 3 used an effective mass 
model to calculate the 15— 15 transition energy Efg^fs 
and Coulomb shifts E^ er for the three CdS/HgS/CdS 
QDQWs. Their findings are also listed in Table III. The 
relative errors between their calculated IS- IS transition 
energies and experimentally measured ones are 8 %, 10 % 
and 11 % for L c i a d — 0.3 nm, 1 nm, 1.5 nm, respectively. 
The ETB theory clearly provides more accurate energy 
gaps than the effective mass model. The predicted gaps 
of the ETB differ slightly from the measured gaps but 
are about 150 meV less than the prediction of the effec- 
tive mass model. Most of these differences are due to 
differences in the single particle energies. The remaining 
part of the difference in energy gaps is due to the differ- 
ence in Coulomb shifts. The Coulomb shifts predicted 
by effective mass theory are 20 meV lower than the ETB 
results. It should also be noted that single-band effective 
theory 3,9 predicts the optically active 1S-1S transition to 
be the ground state transition. Experiment 5 and tight- 
binding theory 11 ' 12 show that the ground state is dark 
and must be the 1S-1P transition as we have calculated 
for the exciton ground state. 

In the tight-bind model, the spin-orbit interaction is 
determined by the parameter Aj = (£i,T| H so \ Zi,\) 
where i = a (anion) or c (cation) and H so is the Hamil- 
tonian of spin-orbit interaction 18 ' 19 . In the bulk, it is 
^ fcn own that the spin-orbit coupling lifts the degeneracy 
coui of the zone-center bulk band states and produces a 4- fold 
■" degenerate state (the light and heavy hole bands) and a 



Table III summarizes our calculated Coulomb shift 
Ecoui, exchange splitting E exc h 1 IS- IS transition energy 
Eis-is an d optical excitonic gap E° pt for the three 
CdS/HgS/CdS QDQWs. E 1S - 1S is the energy of the 
lowest optically active transition because it is a transi- 
tion between the even parity IS electron state and the 
odd parity lSi hole state. E° pt is the energy of low- 
est possible spin-singlet transition which is between the 
lowest electron (even parity IS) and lowest hole (even 
parity lPo). E opt corresponds to the emission peak in 
experiment. In the experiment, the points of maximum 
curvature in the absorption spectra are a good measure 
for the IS- IS transition energies 20 and the experimen- 
tally measured gaps E c g xp3 for the three CdS/HgS/CdS 
QDQWs are listed in Table III. We find that our calcu- 
lations for IS- IS exciton energies are in good agreement 
with experimental E e g xp (the relative errors between the 
experiment and our theory are 0.2 %, 2.7 % and 3.9 % for 
Ldad — 0.3 nm, 1 nm, 1.5 nm, respectively). Also Table 
III shows that our calculated Coulomb shift varies from 
107 meV to 96 meV and the exchange splitting from 1 



2-fold degenerate state (the split-off band) at the zone 
^ center. For example, in bulk CdS, the zone-center split- 
^ tin g between the split-off band and the light and heavy 
hole bands is about 80 meV 10 . Bryant and Jaskolski 12 



have shown that the spin-orbit splittings of single-particle 
states in the QDQWs are much smaller than the bulk 
zone-center ones. At large wavevector, the bulk CdS and 
HgS band structures show little effect of spin-orbit cou- 
pling. Spin-orbit effects and mixing are weak in QDQWs 
because the strongly confined trap states are made from 
bulk states with large wavevector. Coulomb interaction 
doesn't strongly mix different levels, so our results should 
not change much when we include spin-orbit interaction. 
Detailed work about the effect of spin-orbit interaction 
on the fine structure of the excitonic spectrum is still in 
progress and will be reported in a future paper. 

The ETB model we have considered for the QDQW has 
the following specific features: a spherical geometry, a 
particular well thickness and position inside the QDQW, 
no faceting of the surfaces or interfaces, a perfect fee lat- 
tice, specific choices for uncertain material parameters 
(e.g., valence band offset, dielectric constant e ), and wa- 
ter barrier not included. However, the electronic struc- 
ture obtained for the QDQW is determined mainly by 
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the trapping and state symmetry and does not appear to 
depend significantly on a precise choices made specifically 
for QDQW geometry and material parameters. Bryant 
and Jaskolski 12 have test this by studying single-particle 
states in QDQWs with different shapes and different well 
geometries and they have determined the dependence of 
the results on valence band offsets and spin-orbit cou- 
pling. For example, they have found that similar elec- 
tron (hole) states exist for spherical, tetrahedral and cut- 
off tetrahedral QDWSs provided that a similar number 
of cations (anions) occupy a particular region for each 
shape; the splitting between electron (hole) levels has a 
weak dependence on well thickness, well position, valence 
band offset, or as already mentioned, spin-orbit effects. 
The main effect of any of these uncertainties would be 
on the absolute energy position of the ground state tran- 
sition. The position of other transitions, relative to the 
ground state transition, are insensitive to the uncertain- 
ties. Our good agreement with experiment shows that we 
also describe the lowest transition well with our model. 
The computational tests of Bryant and Jaskolski 12 are 
not exhaustive. Other possibilities (e.g., faceting at sur- 
faces or interfaces, dependence on other tight-binding pa- 
rameters, the dielectric screening, or on the effect of the 
water barrier) could be considered. However, further ex- 
perimental characterization of QDQW geometry and ma- 
terial parameters and reduced experimental uncertainty 
are needed to provide tighter constraints for more com- 
plete sensitivity tests. 

In summary, we have studied electron-hole correlations 
in QDQWs by incorporating Coulomb and exchange in- 
teractions into an empirical tight-binding model. Our 
calculated optical excitonic gaps are in good agreement 
with the experiment. Our ETB theory can provide a 
good description for QDQW's nanocrystals. 
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